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Abstract. The dynamical consequences of a bimetric scalar-tensor theory of gravity with a dy- 
namical light speed are investigated in a cosmological setting. The model consists of a minimally- 
coupled self-gravitating scalar field coupled to ordinary matter fields in the standard way through 
the metric: g^ + Bd^dvij)- We show that in a universe with matter that has a radiation-dominated 
equation of state, the model allows solutions with a de Sitter phase that provides sufficient inflation 
■ to solve the horizon and flatness problems. This behaviour is achieved without the addition of a po- 

tential for the scalar field, and is shown to be largely independent of its introduction. We therefore 
have a model that is fundamentally different than the potential-dominated, slowly-rolling scalar 
field of the standard models inflationary cosmology. The speed of gravitational wave propagation 
is predicted to be significantly different from the speed of matter waves and photon propagation in 
the early universe. 

(N 
> 

\£> • 1. Introduction 

(n : 

Recently we have introduced a class of models that could loosely be described as a minimal 
introduction of dynamically evolving propagation speeds for different fields. These models ffl, pL |3| 
use a prior-geometric combination of fields to form a metric that is universally coupled to matter 
fields, thereby not introducing violations of the weak and Einstein equivalence principles. They are 
essentially a dynamical realization of the earlier work |||, |6| by one of the present authors on the 
q-i possibility of a phase transition in the speed of light in the very early universe solving the horizon 
5^ | and flatness problems of cosmology. 

The model considered herein is that introduced in j|] , and uses a scalar field <fi that is minimally 
coupled to a gravitational field described by the metric g^ u , to construct the matter metric: 

X " = + Bd^dvQ. (1) 

It is this metric that is used to construct the matter action and can be said to be the geometry on 
which matter fields propagate. Throughout we will refer to as the gravitational metric, g^ as 
the matter metric, and since (f> is being used to introduce a bimetric structure through (|l]), we shall 
refer to <fi as the biscalar field in order to distinguish it from other scalar fields that may appear in 
the matter model. It is the combination of the gravitational metric and the biscalar field that we 
consider as being the gravitational fields of our theory. 

This type of model has attracted some interest in the field (7|, ||, [|, and similar models that 
introduce the prior-geometric structure in a different way have also appeared 11 1. Since a 
varying light velocity (in 3+1 spacetime) is a fairly generic feature of higher-dimensional models fll"2| , 
13, 14], they are also of interest as effective or toy models that may result from such a theory. 
There is also the more phenomenologicaliy-based 'varying constants' theories fl5|, [l6| , [T7| , [lcf| , and 
investigation of the quantum implications of such models in the very early universe [ ]lS| , [20[ pi]. Lest 
one think that these models are purely of theoretical interest, an analysis of quasar spectra [2^, p3| 
and of the cosmic microwave background [24] have both indicated the possibility of time- variation 



Date: February 7, 2008. 

Key words and phrases. Cosmology, causality, inflation, alternative theories of gravity. 
PACS: 98.80.C; 04.20.G; 04.40.-b. 

1 



2 



M. A. CLAYTON AND J. W. MOFFAT 



of the fine structure constant, and placed constraints on the variation of the speed of light since 



last scattering |g5|, [26|] . 

In this work we will explore the possibility that our model contains cosmological solutions with 
sufficient inflation to solve the horizon and flatness problems. We will begin in Section [2] by 
reviewing the model introduced in || and its reduction to a homogeneous and isotropic spacetime. 
As we shall see, the resulting field equations are quite different than those of standard Einstein 
gravity plus matter or scalar-tensor models, and in Section || we present an (implicitly defined) 
solution for vanishing matter. We find that prior to a time t p t, (massless) matter fields will continue 
to propagate with a velocity c (we will work in a comoving frame of the matter metric g^), whereas 
the speed at which gravitational and biscalar field disturbances propagate, which is determined by 
will become vanishingly small. This 'splitting' of the two light cones causes the universe to 
experience exponential inflation, without requiring the existence of a self-interaction potential for 
the biscalar field. At t p t the restoration of Lorentz symmetry occurs, and the light cones of all 
fields coincide. Assuming that f re h > tpt (the former is the reheating time), we show that we can 
have sufficient inflation to solve the horizon and flatness problems. 

In Section ||] we re-introduce matter with a radiation equation of state, and find an approximate 
solution that, in addition to the symmetry restoration time t p t, has an earlier period (i < t t ) during 
which the light cones of the gravitational and biscalar field are also 'split'. 

One fundamental issue that we will encounter is how to define "the Planck scale" . Combinations 
of quantities such as, for example, £p = y Gh/c 3 , depend on the speed of light c, gravitational 
coupling G and Planck's constant h. Since in our model we have different propagation speeds for 
different fields (which are also dynamically varying) and we can expect that the effective gravita- 
tional coupling to matter also changes with time, we should expect that the scale at which quantum 
effects become important will be different for different fields. It is also normally assumed that h 
is a fixed constant, but clearly if we are allowing fundamental constants to vary with time, then 
it is not unreasonable to consider it to be represented by a field: h(t,x). Although we will not 
completely resolve these issues, we will motivate the choice of the parameter B appearing in (|l]) as 

B » —t\. (2) 
32vr P V ' 

In Section || we show that including matter with a radiation equation of state will not alter 
the de Sitter behaviour significantly, but it will introduce another time scale tt at which the light 
cone of the biscalar field deviates from that of the gravitational field g^ u . We therefore see that 
in general we will have two time scales in the model, each associated with the bifurcation of light 
cones associated with different fields. In Section ^, we end with concluding remarks. 



2. The model 

The model that we introduced in consisted in a self-gravitating scalar field coupled to matter 
through the matter metric (|l]), with action 

S = S giav + Sfj, + iSmj (3) 

where 

S grav = ~jdfi(R[g] + 2A), (4) 

k = 167rG/c 4 , A is the cosmological constant, and we employ a metric with signature (+,—,—,—). 
We will write, for example, d\i = v / — g d x and [i = \f—g for the metric density related to the 
gravitational metric g^, and similar definitions of dfi and jl in terms of the matter metric g^. The 



minimally-coupled scalar field action isQ: 



S* = ~ fdfi [~g^d M <R0 - V{<j>) 



(5) 



where the scalar field <p has been chosen to be dimensionless. The energy-momentum tensor for the 
scalar field that we will use is given by 

1 



+ 



r fiu 



v(</>) 



(6) 



and is the variation of the scalar field action with respect to the gravitational metric: SS^/Sg^ = 

Instead of constructing the matter action §m using the metric g„ v , we use the combination (|l]) 
resulting in the identification of g^ v as the physical metric that provides the arena on which matter 
fields interact. That is, the matter action Sm[iP ] = §M[g, tp 1 ], where ip 1 represents all the matter 
fields in spacetime, is one of the standard forms, and therefore the energy-momentum tensor derived 
from it by 



5Sm 



satisfies the conservation laws 



V, 



0, 



(7) 



(8) 



as a consequence of the matter field equations only Q. It is the matter covariant derivative V M 
that appears here, which is the metric compatible covariant derivative determined by the matter 
metric: V^g^ = 0. 

As described in [01, the gravitational field equations for this model can be written as 



5 2 * 2/z 

and that for the scalar field (written here in terms of matter covariant derivatives) as: 

gT$v$v<l> + KV'[<j>] = Q. 
In the latter, we have defined the biscalar field metric 



yV LV 



K [i 



and 



K = 1 - BeTd, 



In this work we will assume a perfect fluid form for the matter fields: 



nOtj3 



PI 



■ r af3 



(9) 

(10) 

(11) 
(12) 
(13) 



c 2 . 



with g^ u u^u p 

We will now specialize to a cosmological setting, imposing homogeneity and isotropy on spacetime 
and writing the matter metric in comoving form as: 

diag(c 2 ,- J R 2 (t) 7i . 



mil 



(14) 



Note that our earlier publication M contained a typo in equation (5); the rest of the letter is correct as published. 
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with coordinates (t, x l ) and 3-metric jij on the spatial slices of constant time. The matter stress- 
energy tensor (using u° = c): 

f 00 = p, fv = j^jV, (15) 

then leads to the conservation laws (an overdot indicates a derivative with respect to the time 
variable t, and H = R/R is the Hubble function) 

p + 3ff(p+-§) =0. (16) 

Since we are interested in the very early universe, we will assume a radiation equation of state: 

V = \c 2 p, (17) 

which leads to p oc 1/i? 4 . 

It is useful at this point to introduce the following quantities derived from the constant B which 
will appear throughout this work: 

c 2 1 
#1=12^ P B = 2^B> (18) 

the latter comes from H B = ^kc^pb- From the definition (|l^) we have 

and using the relation (Q) the gravitational metric is found to be 

= dmg(Kc 2 ,-R 2 li] ), (20) 
and so p = \f~Kp. From the definition ( |i"l"| ) we find: 

r ' = diag (^'-i^ 7 ^)' (21) 

where we have defined 

P IPB P GPB 



(22) 



In this work, we will find that neither I p nor I p vanish in the solutions of interest, so the metric (|1 
is degenerate only when R = or K = 0. 

The Friedmann equation (i?oo + Kc 2 ^ 13 Rij) is given by 

kc 2 K 1 1 
H 2 + = -Kc 2 A + H 2 B (1 -K) + 2H 2 B KV B + -Kc 4 K 3/2 p. (23a) 

R z 6 6 

The remaining equation from (^) is 

K kc 2 K 1 
2H + 3H 2 - H— + — — = c 2 AK - 3H 2 B (l - K) + 6H 2 B KV B - -nc^^Kp, (23b) 
K R z 6 

which can be shown to be redundant given the scalar field equation (below), the matter conservation 
laws (16) and the Friedmann equation. The scalar field equation fllOD in this gauge is 

I P 4> + 3KHI p <j> + c 2 K 2 V' = 0, (24a) 

which, using ( |l9| ) can also be written as 

I p k - 6K(1 - K)HI P - 2K 2 V B = 0. (24b) 

We have introduced the dimensionless potential: Vb = BV and Vb = V B <j). 



These combined equations suggest the following definitions motivated by the varying constants 
theory fl|, ||]: 



c vc (t) = yfKc, G vc (t) = K 3 / 2 G, (25) 
in which case the Friedmann equation becomes 

H 2 + ^ = \<ZA + ^ + \c 2 KV + f GvcP . (26) 

Aside from the biscalar potential term (which for simple polynomial potentials can be re-written 
as a potential with time varying mass and coupling-constants) we have a form that is identical to 
the Friedmann equation for matter and a scalar field minimally coupled to Einstein gravity. This 
suggests the definition of a Planck scale that is no longer a constant quantity, but varies as, for 
example: 

«■■-<« - M> - js*- (27> 

where pp = c 5 /(hG 2 ). 

We should note that in [^] it is stated that varying constants theories require a violation of 
the strong energy condition in order to solve the horizon and flatness problems. Since we will be 
setting the biscalar field potential to zero, one might ask how it is possible to claim to solve these 
problems with matter that satisfies a radiation equation of state? The answer to this is that the 
Friedmann equation ( |23a| ) does in fact display coupling to sources that appear to violate the strong 
energy condition. Note though that to a material observer, matter fields will obey the strong energy 
condition but couple to the gravitational field in a way that does not. Thus it is the reaction of 
the gravitational field that sees such violations, which is precisely what is necessary for inflation. 

3. EMPTY SPACE, BISCALAR FIELD MODEL 

The simplest possible scenario is one in which there is no matter present, and the biscalar field 
is solely responsible for the dynamics of the very early universe. As we will see in the following 
section, the presence of matter will have a nontrivial effect at earliest times, but not in a way that 
will affect the scenario significantly. 

From we see that both K = and K = 1 are solutions, the first corresponding to an exact 
de Sitter solution and the latter to a model in which the biscalar field plays no active role. It 
is not hard to show that the K = solutions are unstable, with deviations K oc i? 6 , and so we 
end up with a strong relationship between how small K is at the beginning of the universe and 
the amount of inflation achieved. This fixed point corresponds to degenerate g Ml/ and g^ v , and so 
physical solutions will have K > except possibly at a point, which could be considered a physical 
singularity. What we are essentially going to find here are solutions that interpolate between these 
two limits, beginning near K = and making a transition (at roughly t pt ) to a minimally-coupled 
scalar field model. 

Setting p = = p and V B = in ( p4b|) , we find the solution 



K = { l + 2 w) (28) 

where we have defined a scale i? pt to correspond to the end of the inflationary phase, at which 
point we will have K p t = 1/3. This is justified by using (|23| ) and (|24b| ) to derive 

R = RHl(l-K)(l-3K), (29) 

which demonstrates that we have inflation (R > 0) when K < 1/3. The subscript "pt" indicates 
that this is the scale at which the phase transition between an earlier 'broken phase' where the 



6 



M. A. CLAYTON AND J. W. MOFFAT 



propagation speeds for matter and gravitation are very different, and a later 'restored phase' where 
these speeds become approximately equal. 

Using ( p8|) in the Friedmann equation ( |23aD , we find the implicit solution 



^2(2 + y) + 2 



where 



V2(2 + y)-ln v - y c = 6H B (t - t pt ), (30) 

V v 2 \ 2 + y) ~ 2 J 



with y c resulting from fixing the constant of integration so that R = R p t at t = t pt . 
For it! <C Rpt we can expand (|30|) in y to find 



R « i? pt exp[i(y c + ln8 - 2)]e H ^ t ~ t ^ « i? pt e HB (*-H (32) 

and since exp[i(y c + ln8 — 2)] ~ 1.27 we use this approximation right up to t P f In the opposite 
limit (R ^> -Rpt), we find that 



R^R 



pt 



■11/3 



1/3 



-^= + 3\/2#B(t-*pt)J «i2 pt [l + 3^2F B (i-i pt )J , (33) 

where in this case we are making only a slightly larger error since y c /y2 ~ 0.111. 

In the very early universe (t <C i p t) the Friedmann equation has a source that is effectively a 
constant energy density pn, so it is reasonable to assume that ps ~ pp = c 5 /(hG 2 ). At i p t the 
effective energy density has fallen to 2/3 of this value (since K pt = 1/3) and at some scale following 
inflation we expect reheating to take place. It is likely that in order to have reheating we will 
need to consider a non-vanishing biscalar field potential, which for concreteness we will assume 
is a simple quadratic form V(</>) ~ 167r/ _2 2 where I is a length scale characterizing the mass of 
the biscalar field. The dimensionless potential is then Vb = 167rBl~ 2 cj) 2 with 32nBl~ 2 ~ (£p/l) 2 a 
dimensionless parameter. 

Reheating will begin when the biscalar field begins to oscillate about = 0, which is characterized 
by the potential dominating the Hubble damping. It is possible that this could happen prior to t pt , 
in which case the analysis becomes quite complicated, since one would have to consider nonlinear 
contributions to (pi]). If it happens at or after t pt , then we can approximate the scalar field equation 
by 

4> + 3H(j) + 32Tr(c/l) 2 (f>^0, (34) 

and find that reheating begins at a scale R re h where H 2 ch ~ 327T x 4c 2 /(9Z 2 ). From this we see that 
reheating begins at or after t p t if I 2 > 8(327ri?) ~ 8£p. We will leave the details of reheating to 
future work, and assume for simplicity that reheating takes place at t re h = t pt and is instantaneous. 
This means that p Te h, the radiation energy density produced at the end of reheating, should satisfy 
-fr|j(l — K pt ) = |Kc 4 /? re h, which gives yO re h = \pb- With this assumption, we would have the 
reheating temperature not too different from the Planck temperature: T pt ~ T ie h ~ Tp. 
The horizon problem is solved, if the size of the cosmological horizon 

'-w- (35 > 

at some time in the history of the universe is much greater than its value at present ^7|, that 
is: maxf r c (i) 3> r c (to). During inflation (t < t p t) the cosmological horizon is decreasing, and 
it will have a maximum value at earlier times. As it stands our purely classical model can be 
traced back to t = — oo and in that limit r c — > oo. More realistically, the model becomes a good 



7 



approximation to a more fundamental quantum model after a time we will call t qg . We then must 
have r c (t qg ) 3> r c (to), which yields the condition RqHq S> R qg H qg , or 

RqHq R qg H q; 



» l^M . (36) 
RptHpt RptHp t 

Assuming that the universe undergoes AA e- folds in this time, we have R p t/R qg = exp(AA), K qg 
i e -6A^ an( j Hp t /H qg \/2/3. Using standard physics following i pt ~ t re h> this constraint then 
becomes 



AT»ln(^ A /l0 6 ^ eq ) € (32,60). (37) 

The redshift at matter-radiation equality is z cq ~ 4.3 x 10 4 f2/i 2 ~4x 10 4 , and 1 < T sr /Tp < 10 -5 . 
Thus we have the usual statement that M ~ 60 will solve the horizon problem. 

The demonstration that this also solves the flatness problem proceeds in a manner similar to 
that described in ||. Identifying the physically meaningful measure of the flatness of a spacelike 
hypersurface as the 3-curvature measured by a material observer, we have (see also [28|) 



kc 2 

|n-i| = SfH- < 38 > 

We see that this increases without bound as R — > 0, but relating the flatness at t qg and t p t we find 



^ — l|pt RqgH qg 3 

and we see that solving the horizon problem implies that we have also solved the flatness problem, 
as expected. 



Using the definition (19) and the solution (|2q ) we can determine how 'far' the scalar field travels 
between two scales Ri and Rf. 

\A<p\ = 2V3HRf/Ri), (40) 

and so between t qg and t p t we have |A<^>| = 2v / 3A/' ~ 200. Note though, that even if the scalar 
field has a relatively large value at t qg , the effect of the potential is essentially removed by the 
factor involving K that multiplies it. Thus the effect of the potential is delayed until after t pt , 
and the magnitude of the scalar field will have an effect on reheating rather than the dynamics of 
inflation. We therefore see that in this model the biscalar field can in fact 'roll' quite far, but since 
the potential (and therefore the value of <ft) has essentially no effect on the de Sitter behaviour, we 
do not need an extraordinarily flat potential to achieve sufficient inflation. (Because it is essentially 
the nonstandard kinetic terms for the biscalar field that are driving the de Sitter phase, our model 
shares some similarities with pre-big bang string cosmology [p9fl .) This makes us believe that this 
model may have something rather significant to say about the cosmological constant problem, since 
vacuum energy in the very early universe should be irrelevant to inflation in this model. 

On the other hand, we do require that at t qg the value of K is extremely small: K qg « exp(— 6AT), 
which could be viewed as a form of fine-tuning. There are a couple of ways of looking at this problem. 
As an initial conditions problem, we really only need to argue that in some more fundamental theory 
we expect that the universe should enter this classical regime close to a vacuum where the light 
cones are quite different. This is probably not unreasonable, for we have seen in the literature 
dimensionally-reduced models with light cone variability similar to that of our model. Tracing our 
model back to this other vacuum we are forced towards K = 0, and so the issue becomes one of 
how this more fundamental theory reduces to a model similar to ours in the very early universe. 
Alternatively, we can transform the model back to a frame (with time coordinate T) in which the 
gravitational metric g^ v is in comoving form, using 3HbT ~ exp(3ffei). From this we find that 



s 



M. A. CLAYTON AND J. W. MOFFAT 



dt4> ~ T&t4>, and we see that the apparent fine-tuning dt<p is at least partially alleviated by the 
fact that T — > in the very early universe. 

Although we have tentatively identified ps ~ pp, it is not at all clear that this is appropriate. 
The reason is that since we have fields that have a varying propagation speed, it is perhaps more 
reasonable to define the Planck scale as (assuming that h remains constant), for example: 

*-sbw (41> 

where c(t) and G(t) are the effective speed of propagation and gravitational coupling of the field in 
question. We would then use this (or similar definitions for other Planck scale quantities) to argue 
that quantum corrections become important at this scale. 

We know from the field equation (^) that the metric g^ v determines the speed of propagation 
of the biscalar field, which reduces to 

c = ^VKc^ VKc. (42) 

Since the Friedmann equation looks identical to that of a minimally-coupled scalar field, we assume 
that G = G, and from ( |4l| ) we would therefore have pp = K^l 2 pp. Since this vanishes as R — > 
we conclude that we must have a scale above which quantum gravity dominates. Assuming that at 
i qg we have gKC 4 /0p(i qg ) = H B (1 — i£q g ), we find that 

PB = T^PP - < /2 PP - 2~ 5 /V 15 V, (43) 

which would require pb to be incredibly small, if the horizon problem is to be satisfied. The varying 
constants definition (E?|) yields a result that goes in the opposite direction: 

y / i^ qg (l — -Kqg) V^qg 

which would indicate that ps is greatly in excess of the usual Planck scale. Note that this iden- 
tification would also result in quantum corrections becoming important after t qg , which is not 
reasonable. 

On the other hand, because in the gravitational frame the system is an ordinary minimally- 
coupled scalar field, we would perhaps expect that the standard definition of the Planck scale 
would be appropriate-as we have assumed here. That is not to say that this issue is closed. Once 
we include matter into the system, we could legitimately claim that the relevant Planck scale is 
different for matter fields. We expect that a perturbation calculation will not only provide the 
seeds for structure formation, it will also illuminate when quantum corrections will dominate the 
very early universe (and therefore pin down t qg more completely). 

4. The model with non- vanishing radiation 

In the last section we described a solution for which the matter contributions were set to zero. 
In such a limit we could have reverted back to the field equations in the "gravitational frame" , and 
chosen coordinates such that g^ u is of comoving form. We would then find a simple scalar field 
equation: d^cp + 3HdT4> = 0, with solution 8t<P oc 1/-R 3 , using which, the Friedmann equation 
would yield R oc T 1 / 3 . To relate this solution to that which we found in the previous section 
requires making the coordinate transformation: dt = \/l -j- 1 / (3H bT^ dT which puts g^ into 
comoving form. This coordinate transformation leaves the time coordinate essentially untouched 
for T > 1/(3Hb), and for T -C 1/(3Hb) we see that t oc InT, so that whereas T G (0,oo) we have 
t E (-co, oo). If we want to retain primordial matter, this simple solution is no longer available to 



us since the matter stress-energy tensor appears in the scalar field equation (10). We have chosen to 



work directly in a comoving matter frame since the interpretation of the resulting metric is clearer, 
and because the field equations with matter are slightly more manageable. 

We now wish to retain the coupling to matter, choosing a radiation equation of state appropriate 
to the very early universe. From ( |23a| ) we see that if K ~ and K 3 ^ 2 p ~ as R — > 0, then the 
Friedmann equation is still dominated by the biscalar field kinetic terms, and we have the same 
de Sitter behaviour for the scale factor. In fact what we will find is that, at least when the scalar 
field potential may be neglected in the very early universe, the K = fixed point is approached as 
R — > in such a way that \^Kp — > constant. 

Let us begin by parameterizing the radiation energy density p as 

P = p sr (^) 4 , (45) 

where subscript "sr" indicates that the quantity is evaluated at the scale R sr where the biscalar 
and radiation contributions to the Friedmann equation would be equal. Defining the ratio 

*= P -, (46) 

PB 

this implies that we have 

(1 - K sr ) « oKH\ (47) 



which has the small-c solution: 



I P = 1- aKZ, I p = l + -aZ. 



I p k = QK(1 - K)H(1 + iaZ), 
I P Z = -H(l + 3K)Z(1 - oZ). 



(50) 

(51a) 
(51b) 



K ST 1 - 2a + 0{a 2 ). (48) 

Note that the constant a defined in (|46[) , although a very convenient parameter to use in making 
approximations at later times, cannot in any way be considered as a fundamental parameter of the 
model. Since we are interested in retaining the de Sitter phase, it is reasonable to assume that 
t ST ^> £qg , and since we expect that p could possibly be as large as ps at t qg but no larger, then it 
will be redshifted significantly by inflation in order to make a a very small parameter. In turn (fi"8|) 
tells us that t ST > t pt , so that reheating happens well before biscalar-radiation equality. 
It is convenient to define the variable 

*-h*(%Y- < 49 » 

so that we have 



Setting the biscalar field potential to zero, we have the following equations for K and Z: 



The first of these equations tells us that K is a nondecreasing function of time provided that H > 0, 
which is certainly the case at present and from the form of the Friedmann equation ( |23a| ) must also 
be true everywhere in the past (recall that K £ (0, 1)). We therefore see that K should begin near 
K ~ and increase towards K = 1 as the universe expands. Since we expect that R S> -R S r an d 
K ~ 1 at present, Z should be small and decreasing towards zero at present. We also see from ( |51b| ) 
that Z must be decreasing everywhere in the past, and as R — > we must have Z — > 1/c This is 
important, since it shows that near R = we must have K cx R 8 rather than the K oc R e found in 
Section ||. This guarantees that the constant Hg is still the only important term in the Friedmann 
equation as R — > 0. 
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Since Z is an increasing variable, we can write dt = Zdz, and using ([51]) we find 

dK _ QK{\- K){l + \aZ) 
~dZ ~ Z(l + 3K)(l-aZ) ' 
which can be integrated to find (A a is a constant of integration): 

(1 - Kf (1 - aZf 



K 



Z 6 





\ A 1 /—I 


\(- 


) --^yfK 


- ^ Rsr ' 





Using the definition ( f49| ) this becomes 

(1-K) 

and since at R = R sv we have y = 1, we must have 

(1-K sr )(l- l^V^ 
and combining this with (H8|) we obtain 



2 6 ^ > 0. 



> -Tl sr / 



(52) 

(53) 

(54) 

(55) 
(56) 



We shall now seek an asymptotic (in a) solution to (|54|), and, re-summing the first term, we find 



K 



1 + 2,o-( ^ 



(57a) 



which is identical to (pSp up to rescaling. Taking if pt = 1/3 as before, this tells us that R ST /Rpt = 
cr 1 / 6 . Note that this is not a uniform expansion (the higher-order corrections diverge as R — > 0), and 
we perform a series expansion in R/R pt to find (once again re-summing the leading contribution): 



K 



R 

Rsr 



1 

— a 

L2 



t R 



at a value ift 



-v 2 . 



(57b) 
and we 



These two approximations match at the scale Rt ~ \- 
therefore use the form (|57b| ) for t < tt and ( |57a| ) for t t < t < t SI . 

Since we expect that a is a small parameter, we can now determine the evolution of the comoving 
scale factor R throughout the early universe. In fact, the only new feature is that for t < t% we 
have K oc R 8 , which does not alter the de Sitter behaviour of the scale factor since we still have 
K w 1. 

It is also useful to show that for small a we have I p ~ 1 throughout the history of the universe. 
To begin with, we take a time derivative of I p to find (retaining the scalar field potential for the 
sake of completeness) : 



-aKZH[5 -9K + (3- K)aZ] - 3aZK 2 V B - 



(58) 



Near R = we have Z ~ 1/a and we easily see that I p begins near I p ~ 1 and decreases until a 



time t r 



where the right-hand side vanishes (ignoring the potential), at which point we find 

9K m - 5 



aZ n 



(59) 



After this I p will increase again and end up at I p ~ 1 as Z — > 0. Therefore (59) corresponds to 
the minimum value that I p will attain. By combining (|59| ) with ( |54|) it is not difficult to see that 



3.64a 1 / 3 and finally I p (t m ) « 1 - 2.02a 1 / 3 . 



| + Koa 1 /' 3 where K » 0.989. This gives aZ n 
Clearly the correction can be neglected for small enough a. 

The scenario that emerges is the following. Since we have chosen to work in comoving coordinates 
of the matter metric g M1 , , (massless) matter fields propagate with a velocity c throughout the history 
of the universe. For times t < tt we have I p ~ 4/3 and so from (42) we have the speed of biscalar 

Kc. 



wave propagation is (4/3)v Kc, whereas that of gravitational modes contained in g pv is 
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Following this time there is a period (t t < t < t pt ) in which the biscalar and gravitational speeds 
coincide at \^Kc, and later (t pt < t) all three speeds coincide and Lorentz symmetry is completely 
restored. 

It is interesting that it is the radiation pressure (rather than the energy density) that is ultimately 
determining the transition at i t . That is, it is nB^/Kp ~ 1/3 that characterizes the very earliest 
phase of the universe, and, while it is responsible for the second splitting of the light cones, it 
has very little effect on the de Sitter behaviour of the scale factor. Consequently, although the 
matter energy density ( |45[) evaluated at it is pt ~ ^ PbIo and is therefore larger than the Planck 
density, the effective gravitational coupling to matter (from fl23a|) ) is K 3 ^ 2 G. Using (^) we have 
the effective Planck density for matter is pp/K 3 , and evaluated at tt this becomes 4 24 pp/cr 6 , and 
we see that the matter energy density is much less than this. 

5. Conclusions 

We have introduced a bimetric structure in a scalar-tensor gravitational theory coupled to matter, 
resulting in different propagation speeds for gravitational waves and matter waves. This can in 
principle be detected in gravitational wave experiments, such as LIGO, VIRGO and LISA |3C, 31|. 



Within this model we have developed an alternative scenario to the standard inflationary models 
of cosmology, demonstrating the existence of a solution to the horizon and flatness problems without 
resorting to generic slow-roll scalar field potentials. In fact, we have shown that the de Sitter phase 
in the very early universe is largely independent of the introduction of a scalar field potential. Given 
ongoing concern about the cosmological constant problem in cosmology (see, for example [55]), we 
feel that this is a very significant feature of our model. 

The model is characterized by a time t p t before which the light cone of fields in the gravitational 
sector is contracted with respect to matter. (In the presence of matter fields, there will be an 
earlier time tt before which the light cone of the biscalar field bifurcates from that of the rest of the 
gravitational sector.) Prior to t p t the universe expands exponentially, driven by the kinetic terms 
of the scalar field, and is an effect that is attributable entirely to the differing light cones of matter 
and gravity. As time passes through t pt a phase transition occurs, the light cones collapse back 
onto themselves, and the universe is effectively described by a minimally-coupled scalar field. 

Matter is coupled to this extended scalar-tensor gravitational sector in a manner which preserves 
the Einstein equivalence principle. Therefore the matter model lives on a geometry described by 
a metric g^, and we do not expect that our model will be in conflict with any direct equivalence 
principle tests, nor will it violate charge conservation in the manner described in |33[] . The strong 
equivalence principle is violated, and it is the gravitational reaction to the presence of matter that 
is altered. This is precisely what provides the apparent equation of state that violates the strong 
energy condition, despite the fact that the matter and scalar field actions are of the standard forms. 

Due to the stretching of the perturbative scalar field modes in the early universe in the inflation- 
ary period, we expect that a scale invariant microwave background spectrum will be predicted, once 
the initial conditions for the scalar field and the metric are adopted. This issue will be explored in 
a later publication. 
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